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This  report  was  prepared  by  Dynamic  Science  (The  AvSER  Facility), 

A  Division  of  Marshall  Industries,  under  the  terms  of  Contract  DAAJ02- 
69-C-0030. 

The  purpose  of  this  effort  was  to  (1)  document  and  classify  the 
most  hazardous  factors  concerning  airframe  crashworthiness,  (2)  seek 
methods  of  reducing  vertical  decelerations  at  the  floor  level  in 
potentially  survivable  crashes,  and  (3)  seek  design  methods  for  main¬ 
taining  the  “protective  shell"  around  the  occupants  in  an  accident. 

The  contractor  achieved  these  objectives  by  conducting  a  study  of  43 
major  accidents  involving  the  UH-1D/H  aircraft  to  determine  what  types 
of  structural  failure  were  contributing  to  Injuries  in  helicopter 
accidents  and  by  developing,  programming,  and  verifying  a  23-degree- 
of- freedom,  nonlinear  lumped  mass  mathematical  model.  This  model  was 
then  used  in  a  parametric  study  of  the  UH-1D/H  aircraft  to  evaluate 
potential  areas  of  crashworthiness  improvement.  This  report  contains 
a  description  of  the  accident  data  study,  mathematical  model,  para¬ 
metric  study,  full-scale  drop  test,  and  the  results  obtained. 

The  conclusions  and  recommendations  submitted  by  the  contractor 
are  considered  to  be  valid;  however,  the  mathematical  model  developed 
has  definite  limitations,  the  most  critical  limitation  being  tl^at  the 
model  considers  only  vertical  impact  loads  and  therefore  does  not  con¬ 
sider  the  longitudinal  and  lateral  components  that  are  usually  also 
present  in  the  helicopter  crash  environment.  A  second  limitation  i:> 
that  it  would  be  extremely  difficult  to  use  thii.  approach  to  model  and 
analytically  study  the  crashworthiness  of  future  aircraft  designs  with 
any  confidence.  This  is  due  to  the  problems  that  would  be  encountered 
in  attempting  to  predict  the  necessary  weight  data  to  apply  to  the 
lumped  mass  simulation  and  the  spring  constant  data  necessary  to  apply 
to  the  various  springs  that  connect  the  masses  of  the  model. 

It  is  the  intent  of  this  Command  to  expand  this  mathematical  model 
to  include  dynamic  response  to  combined  crash  loading;  i.e.,  crash  loads 
which  possess  vertical,  longitudinal,  and  lateral  components,  thereby 
developing  a  more  realistic  end  useful  analytical^ tool. 

This  report  is  divided  into  two  volumes.  Volume  I  contains  a 
description  of  the  accident  data  study,  mathematical  model,  parametric 
study,  full-scale  drop  test,  and  the  results  obtained.  Volume  II  is  a 
user  manual  for  the  computer  program  developed. 


SUMMARY 


A  mathematical  model  which  may  be  used  to  determine  the 
dynamic  response  of  a  helicopter  airframe  subjected  to  ver¬ 
tical  crash  loading  has  been  developed. 

This  report  is,  in  effect,  a  manual  which  will  facilitate  the 
use  of  the  computer  program  "CRASH".  The  program  was  written 
to  solve  the  equations  and  to  handle  the  nonlinearities  and 
constraints  which  result  from  use  of  the  mathematical  model. 

The  program  was  used  to  evaluate  the  response  of  the  UH-1D/H 
helicopter  to  vertical  impact  loadings .  Recommendations  have 
been  made  which,  when  implemented,  will  reduce  the  forces 
transmitted  to  the  floor  and  transmission  of  the  aircraft. 
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S  INTRODUCTION 

*  H.  i 

This  "User  Manual"  provides  the  details  and  fundamental  con¬ 
cepts  required  to  execute  "Program  CRASH"  in  a  productive 
manner.  Ihe  mathematical  background  for  Program  ^CRASH  is-  con¬ 
tained  in  Vol.  I  of ‘this  report.  ) . 

>  / — *  .  +  \ 

*  „  .  '  V 

Program  CRASH  is  a  computer  simulation  of  a  general  rotary- 
wing  aircraft.  It  represents  an  attempt  to  provide;  the.epgi-'  . 
neer  with  a  practical  tool  for ' investigating  the  dynamic  re¬ 
sponse  of  such  a  structure  when  subjected  to  vertical  crash  €. 
conditions.  ,  '• 

It  .is  very  important  that  the  user  understand  the  model  de¬ 
scription  and  limitations.  As  would  be  .expected,  the  ‘model's 
geometry,  physical  characteristics,  and  structural  properties 
can  be  quite  extensive.-  To  -facilitate  the  user's  understand- 
ing  of  the  program  an  example  problem  is  demonstrated,,  illus¬ 
trating  the  type  of  problem  which  ,  can  be  solved-.  In  this  ex¬ 
ample,  the  generalized  model. is  adjusted  to  rep-resent  the 
Army's  UH-lD/H  aircraft.  The  numerical  input  ds  dis'cussed  in 
detail.  Representative  output  printer  plots,  directly  from 
the  computer,  arte  also  included  to  further  illustrate  the  type 
and  format  available  to  the  user. 

Although  most^ users  will  be  interested  in  the  attenuation  of 
acceleration  levels  at  particular  locations  due  to  the  non¬ 
linear  deformation  of  the  structure,  the  available  output  is 
not  so  limited  in  scope.  Velocities,  absolute  and  relative 
deformations  (both  linear  and  angular),  forces,  bending 
moments,  and  the  concept  of  energy  absorption  are  also  avail¬ 
able  output  features. 

The  output  takes  three  forms:  a  tabulation  of  input  data,  a 
tabulation  of “output  parameters  at  selected  time  increments, 
and  combination  tabulation-plots  of  output  parameters  called 
for  through  the  input.  Dummy  input  coding  sheets  are  also  in-  . 
eluded  to  assist  the  beginning  user  -in  organizing  the  input. ’ 


MATHEMATICAL  MODEL  DESCRIPTION 


GENERAL 


The  mathematical  model  used  in  thi§  study  represents  the  air-  - 
frame  structure  of  a  rotary-wing  type  aircraft.  It  is  a  non¬ 
linear  lumped,  mass,  model  having  23  degrees  of  freedom.  To  con¬ 
struct  the  model,  the  airframe  structure  is  divided  into  four 
vertical  and  three  longitudinal  sections  as  shown  in  Figure  1. 
Individual  masses  are  identified  by  number  in  Table  I.  The 
vertical  section  divisions  are  (1)  transmission,  engine,  and 
rptor  section,  (2)  mass  above  the  floor  section,  (3)  mass  be- 
.  low  the  floor  section,  and„{4)  landing  gear  section.  The  lon¬ 
gitudinal  section  divisions  are  (1)  nose  section,  (2)  central 
section,  and  (3)  tail  section.-  The  14  masses  are  spring  con¬ 
nected  into  the  model  and  are  used  to  represent  various- sec¬ 
tions  of  the  airframe  structure.  Springs  connecting  the  masses 
are  shown  in'  Figure  2  and  identified  in  Table  II.  The  vertical 
section  masses  permit  a  parametric  study  of  the  distribution  of 
the  load- limiting  properties  throughout  the  important  vertical 
sections  of  the  airframe  structure.  The  longitudinal  masses, 
shown  in  Figure  3,  permit  a  study  of  plastic  hinges  and  shear 
failures  at  four  simulated  airframe  locations.  All  masses  may 
not  be  required  to  represent  a  particular  section  of  the  air¬ 
frame  structure  under  study;  however ,.  uhe  model  is  generalized 
to  the  degree  that  single-  as  well  as  multi-engine  aircraft  may 
be  analyzed. 

VERTICAL  SECTION  DESCRIPTION 

The  masses  representing  the  four  vertical  sections,  when  con¬ 
nected  into  the  model,  simulate  load-deflection  characteristics 
They  are  connected  into  the  model  with  16  direct-  and  seven  far 
coupled  springs.  Each  spring  is  combined  in  parallel  with  a 
damper  as  shown  below. 


///////// 


SPRING  — 


/ 


-M  DAMPER' 


Model  Description  (Lumped  Masses 


TABLE  I. 

MODEL  DESCRIPTION  (LUMPED  MASSES) 

Mass  N^t. 

Description  , 

/ 

Ml 

'  M2 

Rotor  Assembly,  Transmission  Assembly, 

/  M3 

and/or  Engine  Assemblies 

M4 

* 

M5 

\ 

M6 

Airframe  Structure  Above  Floor  Level 

M7 

M8 

M9 

Airframe  Structure  Below  Floor  Level 

M10 

and  Floor  Dead  Loads 

Mil 

/ 

/ 

M12 

Ml  3 

Landing  Support  System 

M14 

/ 

.  --  .  --  -  . . . —L 

Two  types  of  spring  damping  are  considered  in  the  model/  in¬ 
ternal  (hysteresis)  damping  and  external  (viscous)  damping. 
Internal  damping  is  introduced  into  the  model  through  /the  load-, 
defection  curve -showu-in— Figure  4 .—  By  assigning  different 
slop^values  to  the  unloading  portion  of  the  generalized  load- 
de f le c t ion~cu r ve/'~_a_hy's ter^s is  cycle  can  be  generated  which 
will  absorb  energy.  The  degre^  to  which  the  load-deflection' 
curves  will  reproduce  the  aircraft  structure  is  dependent ^upon 
the  quality  of  data  available  for  the  particular' aircraft  to  be 
studied,  and  upon  t^ie  user's  ability  to  interpr^tr'available 
data  and  comprehend, the  dynamics  of  thev  deforming  structure. 

External  damping  is ' introduced  by  the  use  of  ,a  constant  applied 
to  the  rate -of -change  of  spring  deformation.  The  numerical 
value  of  this.-constant  may  be  determined  by  an  analytical 


TABLE  II. 

MODEL  SPRING  IDENTIFICATION 

Spring  No. 

Description 

K1 

Direct-Coupled  Load-Deflection 

K2 

Characteristics  of  Rotor/Trans - 

K3 

mission  and/or  Engine  Support 

K4 

System  to  Upper  Fuselage 

K5 

K7 

Direct-Coupled  Load-Deflection 

K16 

Characteristics  of  Airframe 

K17 

Structure  Above  Floor  Leve]  ' 

e 

K8 

K10 

Kll 

K12 

Direct-Coupled  Load-Deflection 

K13 

Characteristics  of  Landing 

K14 

Support  System 

K20 

K21 

/ 

K15 

Direct-Coupled  Load- Deflection 

K18 

Characteristics  of  Airframe  Struc- 

ture  Aboye  Floor  Level  During  Shear 

K19 

Direct-Coupled  Load-Deflection 

K22 

Characteristics  of  Airframe  Struc- 

ture  Below  Floor  Level  During  Shear 

K27 

K28 

Far-Coupled  Load-Deflection 

•* 

K29 

Characteristics  of  Rotor/Transmission 

K30 

and/or  Engine  Support  System  to  Floor 

* 


TABLE  II.  Continued  . 

Spring  No. 

Description 

K31 

Far-Coupled  Load-Deflection 

K32 

Characteristics  of  Airframe  Struc- 

K33 

ture  Belova*  Floor  Level 

T5 

Torsional  Load-Deflection 

T8 

Characteristics  of  Forward 

Fuselage  Section 

T7 

Torsional  Load-Deflection 

T10 

Characteristics  of  Rear 

Fuselage  Section 

estimate,  based  on  theoretical  4ata,  or  by  performing  a  series 
of  computer  runs  using  various  estimated  values  for  these  con¬ 
stants  and  comparing  the  results  to  experimental  data. 

The  extensional  springs  in  the  model  are  classified  into  two 
types,  depending  upon  the  ability  of  the  structure  represented 
to  restrain  a  tensile  rebound.  A  type  1  spring  is  one  which 
.can  restrain  tensile  rebound,  while  a  type  2  spring  cannot. 

Refer  to  Figure  4  and  the  definition  of  SD(I,7). 

* 

LONGITUDINAL  SECTION  DESCRIPTION 

The  masses  simulating  the  three  longitudinal  sections  of  the 
airframe  structure  are  connected  into  the  model  with  four  tor¬ 
sional  springs  and  four  shear-type  springs.  Connection  de¬ 
tails  are  shown  below. 


7 


Figure  3.  Longitudinal  Section  Spring  Connection  Diagram 


Description  of  Load-Deflection  Curve. 
DESCRIPTION 


SD(I,1}' 

i 

SD(I,2j 
SD(I,^) 
SD(I,<i) 
SD (I ,5) 
SD(lj6) 


Slope  of  linear  elastic  portion  of  curve. 
Deflection  which  causes  yielding  to  occur. 
Slope  of  first  plastic  portion. 

Deflection  at  which  plastic  slope  changes. 
Slope  of  second  plastic  portion  of  curve. 
Unloading  slope. 


2  | Spring  type  | 


Type  "1"  follows  curve  1 
Type  "2"  follows  curve  2 


SD(i;7) 

i 

I 

SD  (1,8) 


Proportionality  constant  for  viscous  damping 
(applies  to  entire  curve) . 


The  interface  between  blocks  A  and  B  is  frictionless,  permit¬ 
ting  relative  vertical  displacements.  The  resistance  against 
such  a  vertical  displacement  is  provided  by  a  shear  spring. 

This  simulates  the  possibility  of  a  shear- type  failure  occur¬ 
ring  in  the  fuselage,  the  severity  of  which  is  controlled  by 
the  shear  spring.  The  two  masses  are  connected  to  blocks  A 
and  B  by  frictionless  pins,  thereby  permitting  relative  angles 
to  form  between  the  two  masses.  Resistance  to  such  rotation  is 
offered  by  the  torsion  spring.  This  simulates  the  formation  of 
a  plastic  hinge  in  the  fuselage.  The  generalized  load-deflec¬ 
tion  curve  used  for  all  springs  is  shown  in  Figure  4. 

MODEL  DIMENSIONS 

The  general  coordinates  and  model  dimensions  are  shown  in 
Figures  5  and  6.  These  dimensions  and  physical  properties  are 
consistent  with  the  data  provided  in  the  Input  Description 
section  herein.  If  the  application  of  the  generalized  model  to 
a  specific  problem  does  not  require  all  of  the  lumped  masses, 
the  weight  of  the  omitted  mass  should  be  assigned  a  value  of 
1.0  pound  and  the  dimensions  describing  its  location  omitted. 

If  the  omitted  mass  happens  to  be  one  which'  has  rotational 
properties,  its  mass  moment  of  inertia  should  also  be  assigned 
a  nominal  value  of  one  unit. 

i  i 


PROGRAM  DESCRIPTION  ' 

The  computational  scheme  of  Program  CRASH  is  controlled >by  a 
small  main  program  which  calls  upon  10  subroutines .  Not  only 
does  the  naturf  of  the  problem  permit  its  cUjM^ction  into  in¬ 
put  and  output  phases,  b-ut  even  the  computatipnal  phase  of  the 
problem  is  uniquely  .suited  for  the  building  block  technique 
provided  by  subroutines,  tfhis  type  of  program  structure  also 
lends  itself  to  overlay  techniques  if  the  user  finds  it  neces¬ 
sary.  Nine  of  the  10  subroutines  are  linked  together  by 
numbered  common,  using  approximately  7,600  decimal  locations. 

The  plotting  subroutine  uses  the  formal'  call  sequence  to  trans¬ 
fer  data.  The  source  deck  is  composed  of  827  cards,  requiring 
approximately  24,600  decimal  core  locations.  The  program  as 
listed  in  Appendix  I  is  compatible  with  an  IBM  360-40  computer. 

Data  is  read  into  the  program  in  two  locations.  First,  the 
main  program  is  informed  how  many  individual  ases  of  data  will 
be  processed,  and,  thereafter,  data  is  entered  into  the  program  % 
only  through  subroutine  "READ".-  The  following  section  provides 
a  detailed  description  of  this  subroutine.  The  input  data  can 
be  grouped  into  foui;.  categories: 

INPUT  ■* 

1.  Physical  Properties:  This  involves  the  dimensions, 

mass  distribution,  and  structure  stiffness  properties  t 
of  the  model.  This  group  of  input  comprises  about 
90  percent  of  the  input;  data. 

2 .  9  Crash  Conditions:  The  vertical  velocity  and  pitch 

rates  of  all  masses  must  be  specified  at  impact, 
along  with  the  ground  conditions. 

3.  Output  Requested :  There  are  a  possible  129  printer 
plot  So  which  can  be  requested/'  Of  these  parameters, 
as  many  as  three  may. be  plotted  simultaneously  on  the 
same  plot  for  comparison.  This  is  the  feature  of  the 
program  which  enables  the  engineer  to  analyze  a  large 
amount  of  data.  The  parameters  to  be  plotted  and 
their.  sequence  must  be  specified.  The  details  of 
this  \output  feature  and  the  required  instructions  are 

■  found  in  the  Input  Description  section  herein. 

4.  Numerical  Integration:  The  total  time  duration  of  the 
numeri cal  i ntegrati on ,  T_sx,  and  the  delta  time  incre¬ 
ment,  DT,  must  be  specified.  The  times  for  which  the 
program  stores  data  for  future  plotting  must  also  be 
specified. 
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MODELING  THE  AIRCRAFT 

The  aircraft  is  modeled  by  dissecting  it  into  its  major  s true- t 
tural  or  weight  components.  The  generalized  model/  Figure  1/ 
is  composed  of  a  number  of  lumped  masses,  defined  by  Table  I,  t 
which  are  used  to  represent  the  major  weight  components  of.  the 
aircraft.  The  plastic  hinges  provided  by  the  model  permit  cer¬ 
tain  masses  to  break  and  rotate  with  respect  to  each  other. 

The  location  of  these  plastic  hinges  is  guided  by  the  user  s 
knowledge  of  the  structure's  ability  to  resist  bending  moments, 
or  past  accident  histories  of  postcrash  configurations.  All  oi 
'’the  lumped  masses,  14  in  all,  mawot  be  needed  to  properly  de¬ 
scribe  a  particular  aircraft.  In  those  cases  where  a  mass  is 
omitted  from  the  generalized  model,  its  weight,  CM(I) ,  must  be 
assigned  a  hominal  value  of  1.0. 

The  load  deflection  characteristics  of  the  model  ate . 
seated  *>y  33  nonlinear  springs  as  shown  in  Figure  2  and  Table 
II.  Ibe.  load  deflection  characteristics  are.  assigned  by  the 
curve  illustrated  in  Figure  4. 

Note*  If  a  spring  is  to  life  omitted,  \ the  blank  input  card 
- *  .  for  that  spring  must' be  present  in  the  input  se¬ 
quence.  The  program  checks  the  spring- type  param- 

-  eter,  SD(I,7)  ,  to  see  if  the  spring  is  present. 

Note : 0  If  the  spring-type  SD(i,7)  =  0.0 ,  the- spring  has 
-  been  omitted  from  the  program.'  If.  the- spring-  . 

-  type  SD(1 ,7)  =  1.0,  the  spring  is- unable  tq  hold 
a  tensi’le  load  or  to  elongate.  If  the  spring- 
type  SD (I , 7)  =  2.0,  the  spring  is  capable  of  re- 
sis  ting  a  tensile  load  and  may  elongate. 

The  ability  of  the  simulator  to  duplicate  a  crash  condition  de- 
pendfipon  the  degree  of  realism  that  tile  user  can  incorporate 
into  the  model*  through  the  use  of  realistic  ioad-defi 

nharartpristics.  The  user  is  limited  only  by  his 
ability  to  idealize  desired  nonlinear  load-de flection  proper¬ 
ties  by  the  series- of  straight  lines  provided  by  Figure  4. 

An  additional  feature  of -Program  CRASH  is  to  limit Jj*1®  1 
tion  of  Certain  springs.  This  is  accomplished  by  an  adjust 
££  to  the  elastic  portion  of  the  load-deflection  curve  when 
the  spring  exceeds  the  user-supplied  limiting  deformation 
value  This  requires  the  numerical  procedure  to  be  stopped 
and  re-started  at  a  former  time,  extending  the  program  both  in 
“jUexity  and  computational . time  required  fur  a  rrm  sequence. 
To  help  reduce" the  excessive  time  required,  the  number  of 
startsPrequired  should  be  reduced.  Towarc  this  end,  the  pro 
gram  is  constantly  projecting  50  time  increments-  ahead. 


estimating  with  the  information  on  hand  if  these  conditions 
will  occur.  If  these  projections  indicate  that  the  program 
will  have  to  be  re-started  50  time  increments  in  the  future, 
adjustments  are  initiated  in  an  attempt  to  correct  the  situa¬ 
tion  before  it  happens.  If.  the  original  load-deflection  slope 
is  too  large,  the  projection  will  be  conservative,  thus  de¬ 
creasing  the  slope  for  the  next  iteration.  If  the  slope  is 
too  small,  the  projection  will  indicate  a  future  violation  and 
the  slope  will  be  increased  for  the  next  time  increment.  The 
amount  of  increase  or  decrease  in  this  siop>e  is  provided  by  the 
user.  If  the  user  supplies  a  very  small-  delta  slope  change, 
the  program  may  not- be  able  to  correct  within,  the  available  50 
time  increments,  thus  forcing  a  re-start  condition.  This  con¬ 
dition  will  continue  until  'the  deformation  limitations  are 
achieved.  This  can  develop  into  a  rather  time-consuming  pro¬ 
cess,  especially-"!' f  there /are  more  than  three  springs  being 
limited.  Thper^egihning  uher  is  cautioned  against  excessive  use 
oJ^  this  program/  feature  in  the  early  stages  of  developing  a 
parti cular^mod^l.  Once  an  idea  of  the  model's  general  behavior 
is  obtained,  the  user  should  have  a  better  idea  of  what  re¬ 
strictions  to  impose  in  this  limiting  process.  The  user  is  re¬ 
minded  that  this  technique  does  not  provide  a  "unique"  load- 
deflection  curve  which  restricts  the  spring  deformation  within 
thedesired  range,  since  there  are  many. possible  curves  which 
-Will  satis fy  the  required  condition.  This  technique  determines 
one  of  the  many  possible  answers.  f 

Note;  The  program  is  written  to  handle  up  'to^To^limiting 
springs  simultaneously. 

CRASH  CONDITIONS  j 

♦ 

A  vertical  impact  with  a  positive  or  negative  pitch  rate  is 
permitted.  Each  single-degree-of- freedom  /mass  is  assigned  an 
inpact  velocity  in  ft/sec.  If  the  mass  his  rotational  proper¬ 
ties,  the  mass  is  given  a  vertical  velocity  and  also  a  possible 
•  nonzero  angular  velocity.  Whatever  the  velocity  condition  is 
at  impact,  it  must  represent  a  kinematically  consistent  sys¬ 
tem.  The  example  provided  represents  the  simplest  system 
possible,  namely,  pure  vertical  input,  all  linear  velocities 
being  equal,  and  all  angular  velocities  being  zero. 

GROUND  CONDITIONS  1 

The  ground  is  modeled  by  a  deflection-force  equation  of  fourth 
order  or  less.  For  example f 

Ground  Deflection  =  GD(I,1)  (force)  +  GD(I,2)  (force)^ 

+  GD(I ,3)  (force)3  +  GD(I,4)  (force)4 


/ 


where  I  ='  1,  2,  3  represents  the  three  possible  contact  points 
between  the  fuselage  and  ground,  X21,  *22*  x23»  of  Figure  5. 

The  deflection  of  the  ground  for  the  contact  point  may  be  quite 
high,  since  a  possible  landing  skid  may  produce  high  soil  bear¬ 
ing  stresses .  ' However,  when  the  fuselage  contacts  the  ground, 
these  soil  bearing  stresses  are  drastically  reduced  due  to  the 
large  increase  in  bearing  surface.  The  general  example  of  this 
concept  is  shown  by  the  following  diagram: 


X21'  X23- 


FORCE 


The  parametric  study  (Vol.  I  of  this  report)  was  performed 
assuming  an  infinitely  rigid  ground  condition,  attempting  to 
approximate  the  worst  condition. 


OUTPUT 


.  1  '  , 

Full  digital  output  is  automatic.  This  is  composed  \of  the 
following  informational  each— requested  plot  time:  * 

1.  Generalized  coordinat^  motion ,  both  verticauJL  and- _ _ 

ancpilar,  during  impact,  taken  from  a  zero  reference. 


2.  Vertical  and  angular  velocities  of  each  mass  during 

impact.  \ 

* 

3.  Vertical  and  angular  accelerations  ofxeach  mass  -  . 
during  impact. 

4.  Relative  spring  deformation  of  each  spring  during 
ifipact. 

•1 

5.  The  force  in  each  spring  during  impact. 


There  are  129  parameters*  any  of  which  can  be  requested  to 
appear  in  the  form  of- e  printer  plot  (described  in  tabular  form 
in  the  Output  Description  section  herein) .  Up  to  three  param¬ 
eters  may  be  plotted  simultaneously  for  comparison. 

Note:  If  more  than  one  parameter  is  requested  on  the 

same  plot,  thought  should  be  given  to  the  expected 
magnitude  of  these  parameters  since  the  plotting 
subroutine  automatically  defines  the  scale  factor. 
This  scale  factor  will  be  applied  to  all  the  param¬ 
eters  appearing  on  a  particular  plot.  Large-mag¬ 
nitude  differences  between  the  parameters  destroy 
the  engineering  usefulness  of  this  feature. 
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GENERAL  PLOW  DIAGRAM  OP  PROGRAM  CRASH 


GENERA!  IfLOW  DIAGRAM  OF  PROGRAM  CI^SH  (CONTD) 


GENERAL  DESCRIPTION  OF  SUBROUTINES  *  , 

ITEM  1,  SUBROUTINE  "READ":  All  necessary  data  for  the  analysis 
of  a  crash  configuration  enters  Program  CRASH  through  this  sub¬ 
routine.  Fixed  and  floating  point  data  must  be  formatted  for 
an  15  and  E10.0  specification,  respectively.  The  order  of  the 
read  statements  is  consistent  with  the  Input  Description  (next 
section  herein) .  Coding  sheet  examples  in  proper  format  and  a 
numerical  example  are  provided  in  the  sections  on  Input  Descrip¬ 
tion  and  Example  Problem.  _ , 

ITEM  2,  SUBROUTINE  -"OPT  1 11 :  All  input  data  entered  through  sub¬ 
routine  "READ"  is  printed  out  by  this  subroutine.  For  a  de¬ 
tailed  description  of  its  format  and  a  numerical  example,  see 
the  Output  Description  section  and  illustrations  therein.  This 
subroutine  is  we  of  four  subroutines  ^capable  Of  outputting  in¬ 
formation  to  the  user.  ' 

ITEM  3,  SUBROUTINE  "RUNGA" :  This  subroutine  executes  tbe 
fourth-order  Runge-Kutta  method  for  solving  first-order  ordin¬ 
ary  differential  equations.  The  four  required  functional  eval¬ 
uations  are  obtained  by  calling  subroutine  "FUNCT".  Briefly, 
this  subroutine  accepts  the  values  of  the  dependent  parameters 
and  their  derivatives  at  t  *  t*  and  calculates  the  values  of 
the  dependent  parameters  and  their  derivatives  at  t  *  t*  +  At. 

ITEM  4 ,  SUBROUTINE  ^UNCT" :  This  subroutine  contains  the  equa- 
tions  of  motion  for  the  generalized  model.  The  transformed 
equations  of>  motion  ar^  of  the  form  3  \ 


IX  (I ) 
^ft 


FUNCT (I) 


Vol.  i  of  this  report  prf^sents  a  more  detailed  description  of 
this  set  of  equations.  Ewysically,  this  subroutine  accepts 
the  forces  and  moments  within  each  of  the  springs  of  the  gen¬ 
eralized  model  (Figure  2) sand  calculates  the  acceleration  and 
velocity  of  each  mass.  Tn|«  subroutine  supplies  subroutine 
"RUNGA"  with  required  information. 

\TTEM  5 ,  SUBROUTINE  "LDEF" :  y Subroutine  "LDEF"  is  the  "load- 
deflection"  subroutine. Knowing  the  physical  characteristics 
of  each  spring  of  the  generalized  model,  it  accepts  the  rela¬ 
tive  deformation  and  elongation  rate  of  each  spring  and  calcu¬ 
lates  the  forces  required  to  produce  these  relative  deforma¬ 
tions.  It  tracks  the  load  history  of  each  spring  so  that  hys¬ 
teresis  during  unloading  may  be  achieved.  .  This  subroutine 
supplies  subroutine  "FUNCT"  with  required  information. 


20 


ITEM  6 ,  SUBROUTINE  "DEFORM" :  Sub routine  DEFORM  (deformation) 
contains  the  equations  of  kinematics  of  the  generalized  model. 
Given  the  current  positions  and  velocities  of  the  generalized 
coordinates,  this  subroutine  calculates  the  relative  deforma¬ 
tion  and  elongation  rate  of  each  spring.  This  subroutine  . 
supplies  subroutine  "LDEF"  with  required  information. 

ITEM  7,  SUBROUTINE  "LIMIT”:  Subroutine  "LIMIT"  uses  current 
spring  elongation  rates  and  relative  deformations  to  estimate 
requested  spring  deformations  50  time  increments  in  the  future, 
i.e.,  tfu^ure  =  tcurrent  +  50  At.  If  these  estimated  relative 
spring  deformations  exceed  user-supplied  limits,  adjustments 
in  that  spring's  load- deflection  characteristics  are  made  in  an 
attempt  to  correct  the  projected  excessive  spring  deformation. 
If  these  adjustments  are  not  sufficient  to  correct  the  exces¬ 
sive  spring  deformation,  subroutine.  "LIMIT"  stops  the  numerical 
process,  resets  time  and  the  necessary  initial  conditions,  and 
recycles  the  program.  If  this  occurs,  subroutine  "LIMIT"  in¬ 
forms  the  user  of  the  time  the  program  was  stopped  and  the  re¬ 
start  time.  ,  —  - - - 

ITEM  8,  SUBROUTINE  "STORE":  This  subroutine  stores  all  of  the 
129  possible  output  quantities  at  the  requested  printer splot 
times.  A  list  of  these  quantities  can  be  found  in  the  Output 
Description  section.  _  _ 

ITEM  9,'  SUBROUTINE  "OPT2" :  This  subroutine  first  prints  out 
the  results  previously  set  aside  in  subroutine  "STORE"  with  the 
exception  of  the  ground  displacements  X 2\ ,  *22  >  911(1  x23*  This 
digital  output  is  not  optional.  Second,  this  subroutine  sorts 
out  the  results  as  directed  by  the  user-supplied  "plot  codes" 
and  presents  the  data  in  the  proper  sequence  to  subroutine 
"PLOTT".  This  output  display  is  a  user's  option.  (Further  de¬ 
tails  are  provided  in  the  Output  Description  section.)  Third, 
this  subroutine  sorts  out  the  load-deflection  history  of  those 
springs  which  have  had  their  deformation  restricted  and  pre¬ 
sents  this  information  in  the  proper  sequence  for  digital  out¬ 
put  display.  This  output  display  is  not  a  user’s  option. \ 

/  \ 

ITEM  10,  SUBROUTINE  "PLOTT":  This  subroutine  can  accept  up  to 
three  dependent  one-dimensional  arrays  and  simultaneously  plot 
their  elements  against  an  evenly  incremented,  independent 
array.  In  our  case,  the  independent  array  is  composed  of  the 
plot  times  stipulated  by  the  user  and  the  dependent  arrays  are 
the  desired  calculated  results  as  stipulated  by  the  user- 
supplied  "plot  codes".  Up  to  three  curves  rttay  be  requested 
simultaneously  on  one  plot  (see  Output  Description  section) . 

No  scale  factors  are  required,  as  the  length  of  the  plot  is 
determined  by  the  time  duration  of  the  run  and  the  number  of 
time  increments. 
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INPUT  DESCRIPTION 


The,  detailed  description  of  the  input  which  follows  is  consis-. 
tent  with  the  order  in  which  the  input  is  read  by  the  computer. 
It  .».s_  sjugges  ted— that-  the'reade'r  simultaneously  review  this 
section  with  the  input  example  provided  in  the  Example  Problem 
section.  ‘ 

The  programmer  has  attempted  whenever  possible  to  form  acronyms 
of  the  input  parameters  to  ease  the  confusion.  The  parameter 
name  in  quotes  should  hfelp  th*e  user  in  recognizing  the  com¬ 
puter  notation  acronym7.  All  input  data  is  floating  point  un¬ 
less  stated  otherwise! 


INPUT  TO  MAIN  PROGRAM 


I 


NRUN 


A  fixed  point  number  informing  the-  main  program  of  the 
"number"  of  cases  to  be  "run"  in  an  uninterrupted 
sequence.  This  card  should  be  followed  by  NRUN  number 
of  input  cases,  each  describing  a  particular  crash  or 
aircraft. 


INPUT  TO  SUBROUTINE  "READ" 


XI  (I)  "Initial  X*'  generalized  coordinate  of  the  model  at 
impact.  XI (1) , . . .XI (14)  describe  the  height  labove 
the  ground  in  inches  of  the  X^...Xn,  generalized 
coordinates.  XI (15) ,  XI (16) , . . .XI (20)  represent  the 
•initial  angles  cif  ©j,  ©2'?*-®6  radians  (see  Figure 
5) .  These  parameters  are  not  used  in  the  computation 
scheme  since  we  are  interested  in  relative, coordinate 
motion;  however,  they  are  necessary  for  a  pomplete 
..  description  of  the  model. 

XDI(I)  "Initial  Derivative"  with  respect  to  time^  XDI(l), 


XDI (2) , . . .XDI (14)  represent  the  vertical  impact  veloc¬ 
ity  of  generalized  coordinates  X.,  X2,.../  in  •’* 
inches/  sec,  and  XDI (15),  XDI (167 ,. ..XDI (207  represent 
the  initial  angular  velocity  of  mass  M5JM6,  ...  M^o 
in  radians /sec .  (see  Figure  5).-  For  example,  XDI(2)/  = 
360.0  and  XDI (3)  \*  360.0  would  imply  that  a  possible 
rotor  and  engine  configuration  had  a  360-inch/sec  •>-  • 
velocity  at  impact.  The  user  is  reminded  here  that, 
if  a  pitch  rate  is  desired  at  impact,  t^he  center  of 
gravity  of  the  aircraft  and  desired  pitch  rate  must  be\ 
used  in  a  pre-run  calculation  to  determine  the  XDI (I) 
input.  Non-zero  XDI (15) ,.. .XDI (20)  would  imply  that 
certain  portions  of  the  aircraft  were /bending  with  re¬ 
spect  to  one  another  at  impact,  or  thpt  these  portions 
had  previously  failed  and  acquired  an  angular  velocity. 

'  This  input  must  be  kinematically  consistent. 
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I 


SD(I,J) 


\ 


CL  (I) 

CM  ( if 

% 


Cl  (I) 


D(I) 


DT 


"Spring  Data".  The  J  elements  completely  describe 
the  load-deflection  characteristics  of  spring  I. 

There  are'-a  possible  331 springs,  eact}  requiring  8 
pieces  of  information  for  its.^description.  Refer  to 
Figure  4  for  a  detailed  definition  of  the  J  elements 
in  SD(TyJ).  For  example,  SD(3,1)  =  20,000',  SD(3,2)  = 
0.5,  SD(3,3)  =0.0,  and  SD{3,4)  =  100.0,  indicating 
that  spring;  3  has  an  elastic  load-deflectioh  slope  of 
20,000  pounds/inch  up  to,  a  relative  deformation  of 
0.5  inch,  at  which  time  its  load-deflection  slope  be¬ 
comes  zero.  If  th e-weight  of  mass  3  was  1000\  pounds, 
CM(3)  =  1000. 0,  then  spring  3  would  be  programmed  as 
a  10G  load  limiter.  SD(3,4)  is  set  at  100.0  inches 
as  a  highly  improbable  deformation  level,  thus  de¬ 
fining  the  entire,  load-deflection  curye  for  the  ex¬ 
pected' limits .  Note  that  SD(23,J)  through  and  in¬ 
cluding  SD(26,J)  are  used  to  input  fene  torsional 

spring  data.  .  \ 

»  •  f 

A  "characteristic  length"  in  inches.  Necessary  to 
describe  the  model's  dimensions.  CL(5) ,  CL(6),... 

CL (10)  are  identified  in  Figure  6. 


"Characteristic  Weight"  of^the  I  mass.  CM(1),  CM (2) , 
...CM (14)  represent  the  weight  in  pounds  of  masses 
M(l),  M(2) , . . .M(14)  (refer  to  Figure  1).  If  any  mass 
is  omitted  from  the  model,  its  CM(I)  value  is  assigned 
a  nominal  value  of  one  (1). 

" Characteristic"  moment  of  "inertia"  of  the  I  mass. 

Cl  (5),  Cl  (6)  . ;. Cl  (10)  represent  the  mass  moment  of 

inertia  in  (in.-lb-sec^)  about  their  respective  mass 
centers.  If  a  rotational ^ mass  is  omitted  from  the 
model,  its  CM(I)  and  Cl  (I)  values  are  assigned  a 
nominal  value  of  one  ( 1 ) 

'  /  *' 

Model  dimensions  in  inghes.  D(l)  through  D(4),  D(6) , 
D(8)  through  D(ll) ,  DJ0L6)  ,  D(17)  ,  D(20)  ,  and  D(21)  are 
indicated  by"  lowercase  letters  in  Figure  6. 

"Delta  Time"  increment.  This  is  the  delta  time  incre¬ 
ment  in  seconds,  ud.ed  in  the  numerical  integration 
scheme.  The*  magnitude  of  the  time  increment  depends 
on  the  crash  conditions  and  properties  of  the  model. 

It  is  suggested  that  the  user  use  the  example  of  the 
Example  Problem  section  as  a  guide.  A  value  of  DT  * 
0.0001  sec  (1/10^  millisecond)  was  found  to  be  a  rea¬ 
sonable  compromise  between  time  required  and  conver¬ 
gence  for  this  example. 


TMAX  "Time  Maximum".  The  numerical  integration  process 
v  will  terminate  at  TMAX  seconds. 

■>  • 

PLT  "Plot  Time".  The  first  time,  in  seconds,  that  calcu¬ 

lated  results  will  be. stored  for  output.  . 

DPLT  "Delta  Plot.  Time"..  Results  will  be  stored  every  DPLT 
seconds,  starting  at  PLT.  seconds  and  ending  at  TMAX 
seconds.  For  example,  if  TMAX  =  0.1»  PLT  =  0.001, 
and  DPLT  =  0,002,  the  results  will' be  stored  at  time 
0.-001  second  into  the  impact  and  every  .002  second 
thereafter  until  the  run  is  complete  at  0.1  second. 

o  r 

Note :  Detailed  digital  output  is  automatically 
provided  at  these  specified  times;  there¬ 
fore  ,  the  PLT  and'  DPLT  parameters  must  be 
.specified  even  if.no  printer  plots  are  to 
;be  requested.' 

"Number  of  Plots  Requested".  A  fixed-point  number 
informing  subroutine  "0PT2"  how  many  individual  plots 
have  been  requested-  There  is  no  limit  on  the  size 

of  NPR. 

6 

"Code"  for  "Plotting".  A  fixed-point  array  providing 
instructions  for  plotting.  Each  of  the  129  possible 
output  parameters  is  assigned  a  plotting  code  as  tabu¬ 
lated  in  the  Output  Description  section*  The  dimen¬ 
sion  of  KP (I , J)  is  KP(NPR,3).  Up  to  three  parameters 
may  be  plotted  simultaneously  on  one  plot.  A  brief 
example  would  be: 

.  KP  (3,1)  =  41  % 

KP>3,2)  =42" 

•  0 

KP  (3 , 3)  =  0 

The  third  plot  would  be™  a  comparison  of  the  accelera¬ 
tion  time-histpry  &f  masses  1  and  2.  Notice -that 
KP(3,3)  *  ’6  teilfe  the  computer  that  a  third  possible 
parameter  will  be  omitted..  Additional  examples  are 
presented  in  the  following  section.  Output  Descrip¬ 
tion. 

NIS*  \  "Number  of  Limited  Springs".  NIS  is  a  fixed-point 

\v  number  which  telis  the  program  the  number  of  springs 
Which  have  had  limitations  placed  upon  their  relative 
deformation.  .The  program  is  presently  dimensioned 
•  for  10  limited  springs... 
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NPR 


VKp(I,J) 


IGR  "Is  Ground  Rigid"?  A  fixed-point  number  specifying 

the  ground  conditions 

IGR  =  0  m  Yes  1 

IGR  =  1  =  No 

If  IGR  is  valued  at  1,  the  program  will  expect  the 
necessary  ground  data  to  follow.  If  IGR  =  0,  the 
ground  data  is  assumed  to  be  zero,  thus  producing  an 
infinitely  rigid  ground. 

ISN(I)  .  "Limited  Spring  Number".  A  fixed-point  number  which 
indentifies  the  springs  which  have  had  limitations! 
placed  upon  their  relative  deformations.  The 'dimen¬ 
sion  of  ISN  is  ISN  (NIS).  If  there  are  two  limited 
•springs,  NIS  =  2,  then  ISN(l)  and  ISN (2)  must  be 
specified. 

DFE(I)  "Deflection  Limit".  This  is  the  maximum  allowable 

deformation  in  inches  that  spring  I  will  be  permitted. 

,  There  must  be  NIS  deflection  limits  specified. 

DSD  (I )  "Delta  ^Spring  Data".  If  the  limited  spring  ISN(l) 

exceeds  the  allowable  maximum  deformation  DFE(l),  the 
slope  of  the  elastic  portion  of  its  load-deflection 
curve  SD(ISN,1)  is  changed  to  SD(ISN,1)  +  DSD(l). 
Further  description  of  ISN (I),  DFE(I)  ,  and  DSD.(I)  has 
•been  presented  in  the  Program  Description  section. 

'  _  \ 

GDd^J)  "Ground  Data".  There  are  three  contact  points  be¬ 
tween  the  fuselage  or  landing  gear  and  the  ground 
(see  Figure  5).  GD(l,jr),  GD(2,J),  and  GD(3,J)  repre¬ 
sent  the  respective  J  coefficients  of  a  possible 
fourth-order  polynominal  describing  the  deflection 
load  characteristics  of  contact  points  X2i,  X 22 •  and 
X22.  The  general  form  of  these  equations  is: 

Ground  Deflection  =  GD(I,1)  (force)  +  GD(I,2)  (force)2 

+ ' GD ( I , 3 )  (force)3  +  GD(I,4)  (force)4 

where  the  deflection  is  in  inches  and  the  force  is  in 
pounds . 
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OUTPUT  DESCRIPTION 


Four  subroutines  are  capable  of  providing  output  information 
to  the ''user.  These  four  subroutines  are : 

1.  Subroutine  0PT1 
\  2.  Subroutine  FUNCT 

3.  Subroutine  LIMIT 

4.  Subroutine  0PT2 

i  * 

Detailed  descriptions  of  these  outputs «  follow. 

OUTPUT  FEOM  SUBROUTINE  "OPTl11 :  Output  from  Subroutine  0PT1 
consists  of  a,  tabulation  of  all  input  data.  '  This  tabulation 
is  formatted  to  be  self-descriptive.  This  output  is  not  a_ 
user's  option.;  This  output'  first  reviews  the  plotting  se¬ 
quence  requested  by  the  user.  Notice  in  the  ekample  output 
(Figure  7)  that  the  fifth  plot  has  requested  a  comparison  of  * 
plotting  parameters  8$,  91,  and  92.  This  means  that  the  input 
data  for  the  fifth  plot  was  KP(5,1)  »  88,  KP(5,2)  =  91,  KP(5,3) 
=  92.  Thfs  will  provide  a  printer  plot  compariag  the,  relative 
elongation  of  far-coupled  springs  28,  31,  and  32  (see  Figure  2 
and  Table  III).  The  actual  printer  plot  resulting  from  this 
instruction  is  shown  in  Figure  8.  The  output  of  subroutine 
OPTl  continues  by  discussing  the  limitations  imposed  on  any 
springs,  a  detailed  listing  of  information  about  each  of  the 
14  masses,  and  a  discussion  of  ground  Conditions.  Detailed 
information  concerning  spring  data  follows.  The  notation  is 
defined  at  the  top  of  the  page  with  listings  of  each  of  the  33 
springs.  This  is  illustrated  in  Figures  7  and  9. 

OUTPUT  FROM  SUBROUTINE  "FUN dr" :  This  is  an  error  message  out- 
put  feature  that  informs  the  user  that  he  has  attempted  to 
divide  by  zero.  This  has  occurred  because  either  the  weight 
or  mass  moment  of  inertia  of  one  of  the  lumped  masses  -is  zero. 
It  is  an  input  error,.  as  explained  in  the  Input  Description 
section . 


OUTPUT  FROM  SUBROUTINE  "LIMIT":  If  restrictions  are  placed 
upon  the  allowable  deformation  of  certain  springs ,  the  program 
may  find  it  necessary  to  stop  and  re-start  in  order  to  satisfy 
these  imposed  limitations .  If  this  re-start  cycle  occurs,  the 
user  is  informed  of  the  titne  at  which  the  program  failed  to 
meet  these  imposed  limitations  and  the  time  of  the  re-start. 
This  output  will  occur  for  each  re-start  process.  This  infor¬ 
mation  will  aid  the  user  in  obtaining  an  intuitive  feel  for  the 
severity  of  the  restrictions  he. is  imposing  on  the  problem.  ;- 
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AIRFRAME  CRASHWORTHINESS  STUDY.  MODEL  NO 


Figure  7.  Output  From  Subroutine  "OPT1"  (Tabulation  of  Input  Data) 


Vertical  Velocity  of  Coordinate 


CODE 


TABLE  III.  Continued 
OUTPUT  PARAMETER 


35 

36 

37 

38 

39 

40 

Angular  Velocity  of  Mass  Number 

/ 

5 

6 

7 

8 

/  9 

10 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 
<  52 

53- 

54 

o 

% 

Acceleration  of  Coordinate.  Number 

* 

1 

2 

3 

4  . 

•  5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

'  55 

N  • 

5 

56 

6 

57 

58 

Angular  Acceleration  of  Mass  Number 

7 

8 

59  , 

"\  9 

60 

\10 

61 

~ 

•  1 

62 

• 

2 

63 

3 

64 

' 

4 

65 

5 

67 

7 

68 

70  ' 

Elongation  of  Spring  Number 

8 

\  io 

71 

11 

72 

_____ _ - 

12 

73- 

13 

74 

• 

14  / 

CODE 


TABLE  III.  Continued 


/  OUTPUT  PARAMETER 


Elongation  of  Spring  Number 


iigular  Rotation  of  Torsional  Spring 


Elongation  of  Far-Coupled  Spring  Number 


Force  in  Spring  Number 


OUTPUT^ FROM  SUBROUTINE  "OPT2" :  This  subroutine  first  provides 
a  detailed  digital  tabulation  of  output  parameters.  The  for¬ 
mat  titles  each  of  the  arrays  according  to  the  generalized 
model  (Figures  1  through  5)  and  lists  the  results  at  each  of 
the  requested  plot-times.  This  output  feature  is  not  a  user's 
option. 


If  any  plots  have  been  requested,  the  \user  is  supplied  with  a 
combination  digital  and  printfer-scaled-plot  display  of  the  plot 
code  parameters .  An  example  is  provided  in  Figure  10 .  The 
plot  will  be  titled  according  to  the  generalized  model  and  a 
digital  listing  on  the  left  side  of  the  paper.  The  first 
column  is  the  tine  and  the  following  columns  are  the  digital 
values  of  the  requested  plot  parameters  in  E  format.  Notice 
that  the  title  identifies  the  curve  and  digital  listing  by 
assign'.ng  a  plot;  symbol.  Uiere  is  no  limit  to  the  number  of 
plots  that  can  be  requested. 

\ 

\ 

\ 
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Figure  9.  Output  From  Subroutine  "0PT1"  (Tabulation  of  Input  Data) 
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Figure  10.  Computer  Output  (Transmission,  Engine,  Rotor  at  20  ft/sec 
Vertical  Impact)  /-Vertical  Accelerations  for  UH-1D/H. 
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EXAMPLE  PROBLEM 

The  computer  simulator  just  discussed  was  used  to  investigate 
the  effects  of  vertical  impacts  on  a  UH-1D/H  helicopter.  This 
aircraft,  shown  in  Figure  11,  is  in  widespread  use  by  the  Army, 
Air  Force,  and  Marine  Corps  as  a  tactical  transport. 

Two  major  categories  of  input  data  were  required:  (1)  weight 
data  to  apply  to  the  lumped  mass  simulation  and  (2)  spring 
constant  data  to  apply  to  the  various  springs  connecting  the 
masses. 

Information,  supplied  by  the  helicopter  manufacturer  -was  used 
to  distribute  the  weight  in  the  lumped  mass  model.  This  weight 
distribution  is  shown  in  Table  IV  for  the  empty  aircraft  and 
for  the  aircraft  configured  as  a  trpop  carrier. 

4  i 

The  helicopter  airframe  structure  was  mathematically  repre¬ 
sented  by  the  system  of  lumped  masses  (Figure  12)  superimposed 
on  a  general  side  view  of  the  fuselage.  Comparison  with  Figure 
1  will  illustrate  the  application  of  the  generalized  model  to 
the  UH-lD/H  helicopter., 

Of  the  14  masses  available  in  the  model-,  only  10  were  used  to 
represent  this  aircraft.  Mass  (Ml)  was  not  used,  since  in 
this  aircraft  the  engine,  gearbox,  and  rotor  are  located  in 
close  proximity.  Mass  (M2)  was-  used  to  represent  the  main 
rotor  assembly, /while  the  engine  and  gearbox  were  simulated  by 
-mass'  (M3)..  Mass  (M4)  depicts  t:he  tail  rotor  and  90-degree  ., 
gearbox.  The  upper  portion  of  the  cockpit,  section  was  simu¬ 
lated  by  mass  (M5) ,  while  mass  .(M6)  represents  the  upper  por¬ 
tion  of  the  passenger  compartment  and  the  aft  fuel  cells. 

Mass  (M7)  simulates  the  tail-boom  structure.  The  floor  loads 
'  and  the  structure  below  the  floor  including  the  forward  fuel 
cells  were  represented  by  masses  (M8)  and  (M9) ,  which  were 
located  to  either  side  of  a  possible  break  point  in^  the  fuse¬ 
lage.  Mass  (M10)  was  not  used,  although  the  possibility  of 
using  it  i|n  combination  with  mass  (M7)  to  simulate  crushing  of  - 
the  tail  boom  was  considered.  However,  since  the  tail-boom 
weight  is  isnial,l,  and  sirice  in  severe  crashes  peak  floor  accel¬ 
erations  will  occur  before  tail-boom  crushing,  the  effect  of 
this  mass  would  be  negligible.  Masses  (Mil)  and  (M13)  were 
chosen  to /represent  the  landing  skids.  Those  masses  not  used 
were  assigned  a  weight  of  1  pound. 

Analysis  of  accident  case  histories  involving  the  UH-lD/H  air¬ 
craft  revealed  that  two  possible  locations  existed  for  the 
formation  of  a  plastic  hinge:  one  just  forward  of  the  main 
rotor7 and  another  at  the  juncture  of  the  tail  boom  and  fuselage. 

\ 


TABLE  IV.  UH-1D/H  WEIGHT  DISTRIBUTION 


MASS 

NO. 

A/C  EMPTY 
(lb) 

TROOP  CARRIER 
(lb) 

TOTAL 

(lb) 

Ml 

0 

>  - 

0 

M2 

9  80 

- 

980 

M3 

10  47 

25 

1072 

M4 

60 

j 

60 

M5 

250 

1 

250 

M6  - 

847 

1300 

2147 

M7 

230 

- 

230 

M8 

1306 

1645  , 

2951 

M9 

270  ' 

1100 

1370 

M10 

0 

- 

0 

Mil 

60 

60 

Ml  2  * 

0 

t_- 

0 

Ml  3 

60 

.  ~ 

60 

M14 

0 

- 

0 

TOTALS 

5110 

4070 

9180 

The  lumped  mass  model  was  therefore*  arranged/ so  that  these 
hinges  were  simulated  between  the  masses  at  stations  110.00 
and  243.89  as  shown  in  Figure  12.  ^ 

The  load-deflection  characteristics  of  the  UH-1D/H  airframe 
structure  were  simulated  by  the  combinations  of  springs  shown 
in  Figure  13.  Comparison  with  Figure  2  shows  the  application 
of  the  generalized  model  to  this  specific  case. 

Spring- constant  data  for  the  system  were  estimated  by  analysis 
of  accident  case  histories  obtained  from  the  helicopter  manu¬ 
facturer  and  the  U.  S.  Army  Board  for  Aviation  Accident  Re¬ 
search  (USABAAR) ,  and  by  inspection  of  wrecked  airframe  struc¬ 
tures  at  the  U.  S.  Army  Aeronautical _ Depot  Maintenance  Center 

(ARADMAC) . 

The  main  rotor  and  transmission  i'n~tlie~TJH^lD/H  are  supported 
by  a  sturdy  box  structure  which  ties  directly  to  the  floor. 

To  simulate  this  structure,  spring  (K2)  was  omitted  from  the 
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system  and  far-coupled  spring  (K28)  was  used  to  support  mass 
(M2)  at  the  floor,  so  the  rotor  and  transmission  loads  would 
bypass  the  upper"  fuselage  structure.  The  load-deflection 
characteristics  of  the  transmission  and  rotor  support  system 
(Figure  14)  allow  approximately  1/2  inch  of  elastic  deflec¬ 
tion  at  a  load  of  8,000  pounds  before  the  transmission  sup¬ 
ports  fail. 

The  engine  mass  (M3)  is  supported  by  the  upper  fuselage  section 
mass  (M6)  through  spring  (K3)  whose  load-deflection  character¬ 
istics  are  similar  to  those  of  spring  (K28) ,  as  shown  in  Figure 
14.  The  sudden  decrease  in  the  load-carrying  ability  of 
springs  (K3)  and  (K28)  represents  local  buckling  of  structural 
members . 


The,  displacement  of  the  tail-boom  mass  (M7)  is  controlled  by 
three  springs:  torsional  spring  (T7)  ,  far-coupled  spring 
(K33),'and  shear  spring  (K18) .  The  load-deflection  character-* 
istics  of  spring  (T7) ,  shewn  in  Figure  15 ,  permit  a  2-degree 
rotation  before  a  plastic  hinge  forms.  This  plastic  hinge  can 
then  rotate  to  15  degrees  before  failure  occurs,  and  the 
tail  boom  becomes  incapable  of  resisting  further  rotation.  / 
This  unlimited  rotation . is  controlled  by  spring  (K33)  ,  whose 
characteristics  simulate  the  tail  boom  striking  the  ground 
after  a  predetermined  displacement  of  the  center  of  gravity 
of  mass  (M7) .  Shear  spring  (K18)  is  essentially  rigid  so  that 
no  shear  deformation  occurs  at  the  hinge  point. 

* 

Rotation  of  the  forward  portion  of  fuselage  masses  (MS)  and 
(M8)  about  the  potential  plastic  hinge  at  Station  110  (Figure 
12)  is  controlled  by  torsional  springs  (T5)  and  (T8) .  The 
load-deflection  characteristics  of  these  springs  are  shown  in 
Figure  15. •  Shear  springs  (K15)  and  (K19)  control  shear  defor¬ 
mation  at  the  hinge  point. -  As  with  spring  (K18) ,  these  shear 
springs  allow- no  shear  deformation. 

The  load-deflection  characteristics  of  the  landing  skids  are 
represented  by  two  sets  of  springs :  one  for  the  forward  por¬ 
tion  of  the  skids  and  one  for  the  rear  portion.  Each  set  con¬ 
sists  of  three  springs:  one  f4r-coupled  and  two  direct- coup led 
with  springs  (K8) ,  (Kll) ,  and  (K31)  representing  the  front  por¬ 
tion  of  the  skids,  and  springs  (K13) ,  (K21) ,  and  (K32)  repre¬ 
senting  the  rear  portion.  Each,  of  these  sets  allows  simulation 
of  elastic  deformation,  plastic  deformation,  skid  failure,  and 
ground  contact  of  the  fuselage..  'The  load-deflection  character¬ 
istics  of  these  springs  are  presented  in  Figure  16. 

Consider  the  set  formed  by  springs  (K8) ,  (Kll) ,  and  (K31) . 
Spring  (Kll)  allows  elastic  deformation  up  to  2  inches,  with 
an  applied  load  of  18,000  pounds,  at  which  point  the  spring 
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gure  14.  Load-Deflection  Curve  (Engine 
and  Transmission)  for  UH-lD/H. 


Figure  15.  Load- Deflection  Curves  for 
Torsional  Springs,. 
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Figure  16.  Load-De flection  Characteristics  (Forward 
Landing  Skid)  for  UH-1D/H. 
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becomes  essentially  rigid.  Spring  (K8)  also  allows  elastic 
deformation  of  approximately  2  inches  with  ati  applied  load  of 
18,000  pounds.  This  spring  combination  allows  a  deflection  of 
4  inches  at  a  peak  total  load  of  about  36,000  pounds,  or  about 
2G's  on  the  9,000-pound  aircraft  for  each  gear.  Both  skids 
then  allow  plastic  deformation  at- the  36,000-pound  load  for  an 
additional  7  inches.  At  this  point  the  skids  fail,  the  fuse¬ 
lage  contacts  the  ground,  and  the  influence  of  spring  (K31)  is 
felt.  This  spring  represents  the  interaction  of  the  fuselage 
and  the  impact  surface  and,  therefore,  carries  no  load  until 
the  total  deflection  exceeds  9  inches. 

A  detailed  listing  of  the  input  us'ed  to  simulate  the  UH-1D/H 
helicopter  is  shown  in  Figures  7 , and  9 . 

The  initial  conditions  for  this  example  are  a  vertical' impact 
velocity  of  240  in. /sec  with  all  masses  having  zero  angular 
velocity.  Output  of  the  program  is  also  included.  Figures  7 
and  9  are  self-descriptive  output  of  subroutine  OPT1.  The  com¬ 
plete  digital .output  of  subroutine  OPT2  is  omitted  due  to  its 
length.  Three  examples  of  the  plotting  format  are  provided. 
They  represent  the  requested  plots  1  (Figure  10),  4  (Figure 
17) ,  and  5  (Figure  8)  as  shown  in  the  output  of  Figure  7. 

Notice  that  the  example  has  imposed  restrictions  on  springs  3 
and  28  (Figure  7).  Possible  load-deflection  characteristics . 
which  provide  results  within,  these  restrictions  are  provided 
by  output  (see  Figure  18) I  *  * 

*  ^ 

The  exact  input  required  for  this  example  is  presented  in 

Figure  19  vising  the  master  input  coding  sheets?'  included  in 

Appendix  II. 
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Figure  17.  Computers  Output  (Forward  Floor  Hinges,  Rear  Floor  at  20  ft/sec 
Vertical'  impact}  Vertical  Accelerations  for  UH-lD/H . 
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